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Monolayers of hard dimer colloidal particles consisting of two connected spherical lobes form a de- 
generate crystal (DC) at high area fractions. In this DC phase, the particle lobes occupy triangular 
lattice sites while the connections between lobe pairs are randomly oriented, uniformly populating 
the three crystalline directions of the underlying lattice. In this letter we report on dislocation nu- 
cleation and propagation mechanisms observed in DCs and show that certain particle orientations 
form obstacles blocking dislocation glide. We find that the mean distance between such obstacles 
is only Zexp = 4.6 ± 0.2 lattice constants in experimentally observed DC grains. In large simulated 
DCs with no grain boundaries this average distance is Zaim = 6.18 ± 0.01 lattice constants, and the 
probability of finding larger obstacle separations decays exponentially as p{Z) = 0.37e~^^*''^ . Dis- 
location propagation beyond the obstructions is observed to proceed through dislocation reactions. 
Assuming that such reactions are the only mechanism used to circumvent these obstacles, we con- 
servatively estimate that the energetic cost of separating a single pair of dislocations in an otherwise 
defect-free DC grows linearly with the separation. This is in stark contrast to the logarithmically 
growing separation energy for crystalline monolayers of spheres, hinting that the material properties 
of DCs may be dramatically different than those observed in crystals of spheres. 



The microscopic motion of dislocations plays a crucial 
role in melting 0] and governs numerous macroscopic 
phenomena observed in crystaUine materials, including 
plastic flow, yield, and work hardening 

BllSil- Stud- 
ies of dislocations in colloidal crystals enable direct visu- 
alization of such processes 0, H, IM EH) [l2| , providing 
an illustrative model for addressing fundamental ques- 
tions in statistical physics and materials science. Thus 
far, such studies have focused on crystals of spherical 
particles, whose defect transport mechanisms are well de- 
scribed by existing models [l, 0; IB ■ Advances in colloidal 
particle synthesis techniques have enabled the produc- 
tion of a variety of anisotropic yet monodisperse particles 
isl 14, 1^ 3, 11, 3 3- Dimer particles are a simple, 
fundamental extension of spherical particles and can be 
found in systems ranging from granular piles [23,[2ll to 
paired adatoms in thin film epitaxy [23 |. Furthermore, 
dimers are exceptional since although they are nonsphcr- 
ical, their constituent lobes can nevertheless occupy the 
lattice sites of crystal structures formed by hard spheres. 
The study of ordered phases formed by such particles 
therefore constitutes a natural expansion of the existing 
body of knowledge on crystals of spheres. 

Here we directly examine the mechanisms for disloca- 
tion nucleation and propagation in a crystalline phase 
formed by dense monolayers of colloidal dimer particles. 
This crystalline phase, known as a Degenerate Crystal 
(DC), was first identified in simulations of dimer particles 
23I 24 1 , and is defined by the following two characteristic 
features. First, the individual dimer particle lobes form 
a triangular lattice; and second, the particle orientations 
are disordered, uniformly populating the three crystalline 
directions of the underlying lattice. We find that disloca- 
tion glide in DCs of colloidal dimers is severely limited by 
geometric constraints formed by certain particle orienta- 



tions. This restricted dislocation motion suggests that 
the material properties of DCs may be dramatically dif- 
ferent from those of crystals of spherical particles. 

We observe dislocation motion in DCs comprised of 
hollow, hard dimer particles with spherical lobes of di- 
ameter 1.36 /xm and lobe separation 1.46 /xm. Using sol- 
gel chemistry, the rhodamine-functionalized silica par- 
ticles are templated from dimer-shaped hematite cores 
and are sterically stabilized and suspended in a aqueous 
solution, yielding nearly hard-core interactions. A de- 
tailed description of the particle synthesis is provided in 
the supplementary materials. The synthesis procedure 
produces 95% pure dimer particles with particle polydis- 
persity < 5%. The suspension is pipetted into a sealed 
wedge-shaped cell, and particle area fraction is controlled 
by tilting the cell so that particles sediment into the view- 
ing region, which accommodates a monolayer of particles. 
Before imaging with an inverted microscope, the cell is 
laid flat, allowing local equilibration until the area frac- 
tion is constant over the entire region of interest. The 
insertion procedure for flUing a wedge cell creates small 
air bubbles; when these bubbles move near crystalline 
regions they induce defect formation and transport. 

The observed mechanisms for dislocation nucleation 
and glide are summarized in Fig. [1] Nucleation occurs 
when a single particle, (highlighted in the image with 
a thick black dumbbell) rotates, creating a pair of dis- 
locations (Fig. [T^,b). Glide is observed to occur either 
through a swinging move in which one particle lobe re- 
mains fixed while the other swings into a new crystalline 
position, or via a sliding move in which a particle trans- 
lates along its axis. Swinging shifts the dislocation by 
one crystalline row, while sliding shifts it by two rows. A 
sequence of a sliding move followed by a swinging move 
is shown in Fig.[TJ:,d, where the sliding and swinging par- 



2 



tides have again been highlighted by thick black dumb- 
bells. 
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FIG. 1: (Color online) Before-and-after micrographs illustrat- 
ing observed dislocation nucleation and glide moves. Particle 
lobes have been marked with dots, nearest neighbor bonds are 
indicated by lines, and dislocations consisting of paired five- 
fold and sevenfold coordinated lobes have been highlighted. 
(a,b) One rotating particle (highlighted with a thick black 
dumbbell) nucleates a pair of dislocations. (c,d) A disloca- 
tion glides down by three rows through a combination of one 
sliding move (upper dumbbell) followed by one swinging move 
(lower dumbbell). 



These observed mechanisms resemble similar mecha- 
nisms in crystals of spheres. In such crystals a pair of 
dislocations is created through the displacement of two 
adjacent particles (Fig. [5^). Each dislocation consists of 
one fivefold and one sevenfold coordinated particle and is 
characterized by a Burgers vector. The dislocations glide 
apart through a succession of moves in which each seven- 
fold particle shifts its relative lattice position by moving 
in the direction of the Burgers vector. This process has 
the net effect of producing slip in the region between the 
dislocations, shifting the left side of the crystal upward 
and the right side downward in Fig. [21d,c. 

Guided by our experimental observations, we note that 
dislocation motion in DCs is restricted by the constraints 
of the local particle configuration. A schematic of a 
DC where dimers are represented by black dumbbells is 
shown in Fig. [JJi. Nucleation occurs when a single dimer 
rotates as shown by the arrows on the gray particle. To 
allow for glide, the dimer lobes marked with arrows must 
shift in a manner similar to that shown for spheres. The 
critical difference is that while spheres are free to move 
independently, these lobes are constrained to move in col- 
laboration with their partner lobes. The three types of 
particle orientations relative to the glide direction are il- 




FIG. 2: Schematic of mechanisms for dislocation nucleation 
and glide in monolayers of spheres (a - c) and dimer particles 
(d - f). (a,b) Displacing two spheres (gray) creates a pair of 
dislocations, each containing one fivefold and one sevenfold 
coordinated particle. (b,c) The pair glides apart when the 
sevenfold particles shift parallel to their Burgers vectors (out- 
lined arrows) while all other particles retain their crystalline 
positions, (d) A dislocation pair in a DC is created by ro- 
tating one particle (gray) so that its lobes move similarly to 
the gray spheres in (a). Glide proceeds through the motion of 
the lobes marked with arrows, either by swinging (blue parti- 
cle) or sliding (green particle), (e) Nucleation and glide have 
the net effect of shifting the left side of the crystal upward 
and the right side downward. This slip leaves swinging and 
sliding particles intact, indicating accomodation of the dislo- 
cation glide. The red particles, however, would have to be 
severed by this deformation. Since the colloidal particles in 
our suspensions do not break, the required dislocation motion 
is blocked by such particle orientations, (f) The sequence of 
green and blue particles is a zipper of length Z = 4. This 
zipper length sets the maximum separation attainable using 
glide. 



lustrated by the green, blue and red particles in Fig. 
Green particles can shift both lobes in the desired direc- 
tion by sliding, while blue particles can shift one lobe by 
swinging. Particles like these therefore permit dislocation 
glide, in concurrence with the experimental observations 
of swinging and sliding moves (Fig. [lt,d). The red par- 
ticles, however, would need to be broken by the relative 
shifting of the crystal rows during the slip caused by glide 
(Fig. [21i,e). Since the particles in our suspensions do not 
break, the orientation of such particles blocks the motion 
of dislocations. Consequently, only sequences of consecu- 
tive swinging and sliding particles allow continuous glide. 
Since their glide-permitting motion is reminiscent of rear- 
rangements in random square-triangle tilings, we define 
such sequences as 'zippers' [25|. In Fig. [^f, we highhght 



a single zipper. The zipper lobes, whose motion enables 
glide, are marked with a 'z', and the zipper length, Z, is 
defined as one plus the number of zipper lobes. 

In crystals of spheres, dislocations can glide arbitrarily 
far apart, but in DCs the zipper length defines their max- 
imum glide separation. Consequently, dislocation mobil- 
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ity in DCs is determined by the distribution of zipper 
lengths. The ensemble of all particle orientations that 
allow glide for dislocations produced by a clockwise rota- 
tion of one particle is shown in the inset of Fig. [31 Naively, 
one might expect to find a zipper of given length with 
probability p{Z) oc (2/3)'^ since 2/3 of the particle orien- 
tations correspond to swinging or sliding moves. While 
this crude approximation accurately predicts that long 
zippers rarely occur, it overlooks important correlations 
between neighboring particle orientations. 




FIG. 3: (Color online) Probability distribution of zipper 
lengths in both experimental (empty squares) and simulated 
(solid triangles) DCs. Counting statistics determine the error 
bars, which for the simulations are smaller than the symbols. 
The average zipper length measured from experimental DC 
grains is Z^^p = 4.6 ± 0.2. Zippers in simulated crystals with 
10'' lattice sites are only slightly longer: Zsim = 6.18 ± 0.01. 
The dotted line is the best fit exponential for Z > 6, having 
the form p{Z) = 0.37e~^/*'*. The inset shows the ensemble 
of glide-permitting particle orientations given nucleation via 
clockwise rotation of the gray particle. Particle configura- 
tions including a subset of these orientations enable glide via 
swinging (blue/dark gray) or sliding (green/light gray). 



Precisely accounting for these correlations is theoreti- 
cally challenging; instead we directly measure p{Z) from 
experimentally observed DCs. Zippers in DC grains are 
measured by randomly selecting a particle and counting 
the number of zipper lobes extending from it. We find 
that, on average, zippers are Z^xp = 4.6 ±0.2 lattice con- 
stants long. [31| The mean diameter of the observed DC 
grains is 10 ± 1 lattice constants. Clearly, p{Z) could 
be affected by this length scale. To determine the zip- 
per length distribution independent of grain size, we pre- 
pare ensembles of large DCs using numerical Monte Carlo 
moves similar to those described in 2j| . The simulations 
generate crystals with 10^ lattice sites, but the mean zip- 
per length is still only Zs^m = 6.18 ± 0.01 (Fig. E]). The 
tail of the simulated distribution is well characterized by 
the curve p{Z) = 0.37e~'^/*'^, in agreement with predic- 



tions of exponentially decaying orientation correlations 
for dimers on a triangular lattice 2^, 27 1. 

While zippers in DCs are on average only several lat- 
tice constants long, shearing or melting processes typ- 
ically require the transport of dislocations over much 
larger distances. Our experimental observations reveal 
a mechanism for surpassing the zipper length limit via 
dislocation reactions. In such reactions two dislocations 
may combine or one may split apart so long as the sum 
of the Burgers vectors is conserved. These reactions al- 
low dislocations to hop onto nearby zippers intersecting 
the glide path but oriented along a different crystalline 
axis. An experimentally observed dislocation reaction is 
illustrated in Fig. H) In this sequence, a dislocation glid- 
ing down from the upper right approaches the end of its 
zipper. The defect undergoes a reaction and splits into 
two new dislocations. One dislocation's Burgers vector 
is aligned with the horizontal crystalline axis and glides 
to the left along an available zipper, while the second 
dislocation has moved to the lower right through a set of 
moves that are slightly complicated by the presence of a 
nearby grain boundary (Fig. |4|d). 




FIG. 4: (Color online) An observed dislocation reaction al- 
lowing a dislocation to hop from one zipper to another. Only 
the relevant defects have been highlighted, and their Burgers 
vectors are indicated by outlined arrows, (a) A dislocation 
gliding down from the upper right is one lobe from the end 
of its zipper, (b) The dislocation has reacted and proceeds 
by gliding down another zipper extending horizontally to the 
left. A second dislocation, visible to the lower right of the 
reaction site, was emitted to conserve total Burgers vector. 

Such reactions could enable dislocations to separate 
over arbitrarily large distances along a pathway of in- 
tersecting zippers. Nevertheless, the existence of such a 
pathway does not guarantee that dislocations in DCs are 
as mobile as those in crystals of spheres. To elucidate 
the difference between the dislocation transport energet- 
ics in the two systems, we compare the cost of separating 
a single pair of dislocations over N lattice constants along 
the direction parallel to their Burgers vectors in an other- 
wise defect-free crystal. In crystals of spherical particles 
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this energy increases as = 2^^°_,^) ln(iV), where a is 
the lattice constant, u is the 2-D shear modulus and v 
is the poisson ratio [3, |2^ . For dislocations separating 
via intersecting zippers in DCs, each dislocation reaction 
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between zipper segments requires both a core energy to 
create new defects and a se par ation energy as one de- 
fect glides along the zipper [28|. The energetic cost of 
separating two dislocations by Na along their Burgers 
direction using the shortest pathway of connected zip- 
per segments with length Zqq increases linearly with N: 

Edc = 2^F(T-^HZo) ( ^^zif° ) • A detailed calculation 
of this separation energy is provided in the supplemen- 
tary matierials. For crooked or fractal-like pathways this 
energy may increase as an even higher power of N . While 
dislocation reactions in which two defects merge can re- 
lease energy, the presence of additional dislocations in 
the crystal does not guarantee that these would combine 
with defects produced at zipper junctions, as would be re- 
quired to lower the energetic cost of separation. Further- 
more, even though vacancy-mediated climb could be used 
to bypass certain obstacles, vacancy transport in DCs can 
only occur via sliding or swinging particle moves, and 
consequently dislocation climb is also restricted in DCs. 



We claim that the material properties of DCs will be 
strikingly different from those of crystalline spheres. DCs 
will be more resistant to plastic flow since dislocation 
glide cannot proceed along a straight line, as is required 
for slip. Furthermore, if the separation energy does grow 
linearly with N, we speculate that this will have impor- 
tant implications for melting. In the KTHNY theory 
of 2-D melting, the crystal to hexaticphase transition 
requires dislocation pair unbinding [H, (ij . The competi- 
tion between the energetic cost of dislocation separation 
and the entropic contribution to the free energy, both 
of which increase as In(iV) for crystals of spheres, de- 
termines a unique melting temperature. If in DCs the 
separation energy increases as N, dislocation unbinding 
may no longer be feasible at any finite temperature. This 
suggests that melting in DCs may occur via additional 
mechanisms. Furthermore, the observed geometric re- 
strictions in DCs may also apply to other dimer systems, 
such as lipids with dimer-like head groups (29j and gran- 
ular dimers 13, 2^- For example, these restrictions help 
explain why avalanches in 2-D piles of dimer beads oc- 
cur at relatively high critical angles and require tumbling 
rather than collective slip 2l|. Additional comparative 
studies between crystals of spheres and DCs should fur- 
ther elucidate how the seemingly benign act of pairing 
particles into dimers introduces constraints that dramat- 
ically alter the material properties of the crystal. 
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SUPPLEMENTARY MATERIALS 

Particle synthesis 

The colloidal dimer particles are templated from sac- 
rificial ck-Fe203 (hematite) core particles prepared via 
aging of a condensed Fe(0II)3 gel as described in [sO]- In 
a typical synthesis, 50 mL of FeCla (2.0 M) is mixed with 
45 mL of NaOH (6.0 M) and 5 mL of the shape-modifier, 
Na2S04 (0.6 M), before being placed in an oven set at 
100°C for 8 days. A layer of rhodamine isothiocyanate- 
modified silica is then coated on the Q;-Fe203 cores us- 
ing a base catalyzed sol-gel reaction under sonication 
at 30°C for 4 hours. A reaction mixture comprised of 
0.4% hematite powder, 0.25% rhodamine dye solution, 
7.1 M deionized water, 0.92 M ammonia, and 17.4 mM 
tetraorthoethylsilicate added batch-wise to isopropyl al- 
cohol medium produces fluorescent silica shells approx- 
imately 65 nm thick. The shells are stabilized with a 
non-ionic surfactant, polyvinylpyrrolidone (30000 molec- 
ular weight), and finally the hematite cores are selectively 
etched by dissolution in 18% hydrochloric acid solution at 
room temperature. Before use, the colloidal suspension is 
titrated to pH 7, washed via repeated centrifugation and 
decanting, and redispersed in deionized water. This syn- 
thesis protocal routinely produces 95% pure dimer par- 
ticles. 



Dislocation separation energy calculation 

In crystals of spheres, nucleated pairs of dislocations 
with equal and opposite Burgers vectors ±5 separate by 
gliding apart along a straight line parallel to their Burg- 
ers vectors (Fig. 31^) ■ This separation allows for re- 
laxation of externally applied shear stresses, and also en- 
ables the unbinding of dislocation pairs, which is a crucial 
component of the KTHNY model of 2-D melting. Using 
continuum models of dislocation interactions, the ener- 
getic cost of such a separation over N lattice constants 
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is Eg = 2t!{i-v) ^"-(-^)i where a is the lattice constant, 
is the 2-D shear modulus and v is the poisson ratio [2^. 
In addition to this continuum interaction, there is also a 
core energy required to create the dislocations. 

Similarly, the cost for separating a pair of dislocations 
with Burgers vectors ±6 in an otherwise defect-free DC 
can be estimated for any particular intersecting zipper 
pathway. For example, consider the pathway schemati- 
cally depicted in Fig. In this scenario, two dislo- 
cations nucleate and initially glide apart along a straight 
line parallel to their Burgers vector, mimicking the sit- 
uation observed in crystals of spheres. However, after 
separating by a characteristic distance Z^a, equal to the 
average zipper length, the defects each undergo disloca- 
tion reactions, turning to travel along a new zipper seg- 
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FIG. 5: Pathways for dislocation pair separation in crystals 
of spheres (a) and DCs (b,c). Pair nucleation occurs at the 
location marked with a star in each pathway, (a) Disloca- 
tion pair separation is achieved through nucleation and glide 
along a straight line parallel to the direction of the defects' 
burgers vectors ±6. This process results in a pair separation 
of Na. (b) In an otherwise defect-free DC, an identical dis- 
location pair ±6 separates by Zqu, the average zipper length, 
before undergoing dislocation reactions and traveling along a 
sequence of tilted path segments of length Zoa. This pathway 
also results in a separation of Na along the burgers vector di- 
rection, (c) The shortest pathway with segment length Zoa 
yields a separation of Na along a direction rotated by 7r/6 
from the axis of the burgers vectors ±6'. 



ment of length Z^a extending along a different crystalline 
direction. At each reaction site, an extra dislocation is 
produced, which is assumed to remain stationary. 

The energetic cost of separating the dislocations along 

2 

the initial straight segment in Fig. ^Sh is 2J(i-iy) i^(-^o)- 
Again, there is also an additional core energy required 
to create a new dislocation at the junction. Each addi- 
tional segment requires more energy, and includes con- 
tributions from interactions between all the dislocations 
along the pathway. The magnitude of the force between 
dislocations decreases as the inverse of their separation. 
Consequently, the largest force on a dislocation moving 
on a given segment will come from the dislocation near- 
est to it. The energy required for such a separation is 
Ereac = 2 2/(1- v) ^"(^o)- Adding up these energy contri- 
butions from each reaction, as well as the cost of the first 
straight segment, the total energetic cost for the pathway 
depicted in Fig. ^St* is then 



Enc{Z,,N) = -^ln(Zo)+ '^^~^°^ ^.eac 
Z7r(I — v) oZq 

This separation energy grows linearly with final separa- 
tion N , in contrast to the logarithmically growing energy 
for dislocation separation in crystals of spheres. We also 
note that this is a conservative estimate that excludes 
the additional core energy costs at each junction. In the 
limit that the average zipper length Zq approaches the 
separation distance iV, the energetic cost obtained for 
crystals of spheres is recovered. 

The above energetic cost estimation is only strictly 
valid for the particular scenario depicted in Fig. ^SJd. 
However, it represents a conservative energetic cost es- 
timate for any pathway spanning N lattice constants 
along the initial Burgers vector direction, since any other 
option would require additional segments and disloca- 
tion reactions. If the defects are allowed to separate 
by Na along any direction, then the shortest possi- 
ble pathway is the — 7r/3,+7r/3 sequence shown in Fig. 
33};. The energetic cost for separating along this path 

is l^lj^HZo) ( '^^f" ), which stiU grows linearly 
with N . More crooked pathways with many more seg- 
ments could result in an energetic cost of separation that 
grows as a higher power of N . 

We note that in a real system, it is unlikely to find 
one isolated dislocation pair separating in an otherwise 
perfect crystal. The presence of a thermal bath of dislo- 
cations may enable other less expensive mechanisms for 
pair separation. However, since all dislocation motion 
is geometrically constrained, the presence of additional 
dislocations in the crystal does not guarantee that these 
would be able to combine in merging dislocation reac- 
tions, as would be needed to lower the energy of sepa- 
ration. Future studies of the mechanisms of dislocation 
unbinding should help elucidate which additional mech- 
anisms can make such separation energetically feasible. 
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